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We describe a new formulation of the functional renormalization group (RG) for interacting 
fermions within a Wilsonian momentum-shell approach. We show that the Luttinger-Ward func- 
tional is invariant under the RG transformation, and derive the infinite hierarchy of fiow equations 
satisfied by the two-particle-irreducible (2PI) vertices. In the one-loop approximation, this hierar- 
chy reduces to two equations that determine the self-energy and the 2PI two-particle vertex 
Susceptibilities are calculated from the Bethe-Salpeter equation that relates them to ^^-^K While 
the one-loop approximation breaks down at low energy in one-dimensional systems (for reasons that 
we discuss), it reproduces the exact results both in the normal and ordered phases in single-channel 
(i.e. mean-field) theories, as shown on the example of BCS theory. The possibility to continue 
the RG fiow into broken-symmetry phases is an essential feature of the 2PI RG scheme and is due 
to the fact that the 2PI two-particle vertex, contrary to its IPI counterpart, is not singular at a 
phase transition. Moreover, the normal phase RG equations can be directly used to derive the 
Ginzburg-Landau expansion of the thermodynamic potential near a phase transition. We discuss 
the implementation of the 2PI RG scheme to interacting fermion systems beyond the examples 
(one-dimensional systems and BCS superconductors) considered in this paper. 
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I. INTRODUCTION 

The two-particle-irreducible (2PI) formalismiiSiiiii 
was first introduced in condensed-matter physics as a 
means to systematically set up self-consistent approxi- 
mations that satisfy conservation laws. It can be cast 
in a variational framework where the thermodynamic 
potential F is expressed as a functional of the single- 
particle Green functioni^ F is essentially determined by 
the Luttinger-Ward (LW) functional $^ given by the 
sum of the 2PI Feynman diagrams. <& is also the gener- 
ating functional of the self-energy and higher-order 2PI 
vertices. The so-called ^-derivable approximations are 
based on truncations of the diagrammatic expansion of 
<& that retain only a finite number or a sub-series of dia- 
grams. They are thermodynamically consistent and sat- 
isfy conservation lawsiii While most applications of the 
2PI formalism to interacting fermion systems have been 
limited to the Hartree-Fock level, recent developments, 
motivated by the physics of high-temperature supercon- 
ductors, have incorporated exchange of spin fluctuations 
within the Hubbard model^SiSJ*' Following general ideas 
put forward by Wetterichjii the aim of this paper is to 
discuss a 2PI formulation of the renormalization group 
(RG) approach to interacting fermion systems. 

The RG has proven a powerful approach for study- 
ing low-dimensional fermion systems, providing a sys- 
tematic and unbiased method to study c ompet ing insta- 
bilities in the weak-coupling limit (Refs. Il2l4(l) . One of 
its main successes has been to explain how unconven- 
tional superconductivity can occur at low temperature 
in systems like organic conductorsiSiiiiiiifiiifi or high- 
temperature superconductor92SiSii2SiS2i24 where the dom- 



inant electron-electron interactions are expected to be re- 
pulsive and favor antiferromagnetism. Although the RG 
can be implemented in different ways, most approaches 
rely on the so-called one-particle-irreducible (IPI) RG 
scheme or variants thereofi^SiSSiil The IPI RG scheme is 
based on an exact RG equation for the generating func- 
tional of IPI vertices. The existence of a Fermi surface 
implies that the interaction amplitudes strongly depend 
on the momenta of the interacting fermions, which leads 
to functional RG equations for the IPI vertices. For this 
reason, most RG calculations in fermion systems have 
been limited to one-loop order and are thus restricted 
to the weak-coupling limit. Another severe limitation of 
the method comes from the difficulty to access broken- 
symmetry phases. Long-range order is signaled by a di- 
verging flow of certain IPI vertices and susceptibilities 
at a critical energy or temperature scale, below which 
the RG equations cannot be continued. It has been pro- 
posed to circumvent this difficulty by introducing an in- 
flnitesimally small symmetry-breaking component in the 
initial condition of the RG equations, or combining the 
RG technique with a mean-fleld approximation at low 
energy42ii^ Alternatively, one can introduce a partial 
bosonization of the action yielding a description in terms 
of both fermionic and collective Hubbard-Stratonovich 
£gj^gjiij32,38^ It will be interesting to further explore the 
applicability of these new RG schemes in various models. 

The main purpose of this paper is to show that broken- 
symmetry phases can be studied with a RG scheme where 
the basic quantities are the 2PI vertices rather than their 
IPI counterparts. To understand this issue in simple 
terms, let us consider the IPI two-particle vertex in the 
spin singlet particle-particle channel obtained by sum- 
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ming the bubble diagrams (random-phase approximation 
(RPA)), 



7 ^ 



70 



1 + ^70 
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where I ^ ln(Ao/A) comes from the non- interacting 
particle-particle propagator. Aq is a high-energy cutoff 
(e.g. the bandwidth), A an infrared cutoff which can be 
identified with the temperature, and 70 the bare (dimen- 
sionless) interaction. In the RG framework, 7 can be 
interpreted either as the IPI vertex of a theory with in- 
frared cutoff A or as the effective interaction of the Wilso- 
nian action with ultraviolet cutoff .\. " Eq. Q becomes 
a differential equation (see e.g. Ref. |4^, 



^7 
~dl 



(2) 



For an attractive interaction 70 < 0, a divergence occurs 
at the energy (or temperature) scale A^ = Aoexp(l/7o), 
signaling the formation of Cooper pairs and the appear- 
ance of long-range superconducting order. In the BCS 
(mean-field) theory, this divergence is cured below Ac by 
the presence of a finite gap in the fermion excitation spec- 
trum. However, in the IPI RG scheme, it prevents the 
fiow to be straightforwardlySfi continued into the broken- 
symmetry phase. 

Let us now reconsider Eq. Q from a different point 
of view. This equation can be seen as a Bethe-Salpeter 
equation in the particle-particle channel with the bare 
interaction 70 as the 2PI vertex (the reason for this 
notation will become clear in Sec. El- Only the non- 
interacting particle-particle propagator ~ I ~ ln(A/Ao) 
is scale dependent in Eq. (QJ, while the 2PI vertex 
is invariant under the RG transformation. 



dl 



0. 



(3) 



Within the BCS theory, the 2PI vertex <i>(^) is therefore 
not sensitive to the transition into the superconducting 
phase. The appearance of long-range order is expected 
to induce an anomalous (i.e. symmetry violating) self- 
energy below Ac. In more complicated cases, where sev- 
eral types of fiuctuations may compete together, we can- 
not exclude the appearance of singularities in the 2PI 
vertices. For instance, in a conductor close to an an- 
tiferromagnetic instability, singularities in the particle- 
particle channel (i.e. in the particle-particle component 
of $(2)) may be induced by nearly divergent spin fiuc- 
tuations. We shall discuss this point in the concluding 
section and explain how these singularities can be con- 
trolled by a proper parameterization of the vertex. All 
these considerations suggest to use a RG scheme where 
the basic objects are the 2PI vertices. 

The outline of the paper is as follows. In Sec. Ill Al 
we briefiy recapitulate the 2PI formaHsm. By means of a 
Legendre transformation, we express the thermodynamic 
potential (grand potential) as a functional T[G] of the 



Green function; we then define the LW functional $[G] 
and the 2PI vertices In Sec. Ill BL we describe the 

RG procedure. We derive the differential equation satis- 
fied by the thermodynamic potential and show that the 
LW functional is invariant under the RG transformation. 
We then deduce the infinite hierarchy of fiow equations 
satisfied by the 2PI vertices. A one- loop approximation is 
then introduced by truncating this hierarchy and approx- 
imating the 2PI three-particle vertex in terms of the 2PI 
two-particle vertex ^'^^^ fSec. Ill C")) . We discuss the con- 
nection between the one-loop equations and their coun- 
terparts in the IPI RG scheme. In Sees. Ill Dl and HTEl we 
give explicit expressions of the one-loop fiow equations in 
the normal phase and discuss the calculation of response 
functions. Some of the general results of Sees. Ill AIII El 
have been previously obtained by Wetterichfii sometimes 
in a slightly different formulation, with a few important 
differences that we shall mention. In Sec. Ill Fl we show 
how the normal phase RG equations can be used to de- 
rived the Ginzburg-Landau expansion of the thermody- 
namic potential in the vicinity of a phase transition. This 
is achieved by considering the Legendre transform F[E] 
of the LW functional $[G], which allows one to express 
the thermodynamic potential as a functional of the self- 
energy. The latter is then split into a normal part and an 
anomalous (i.e. symmetry violating) part which is to be 
determined by minimizing the thermodynamic potential. 
The last two sections are devoted to the application of the 
2PI RG formaHsm to two different models. In Sec. IIIII 
we consider a three-dimensional fermion system with an 
attractive interaction. From the one-loop approximation 
restricted to the particle-particle channel, we rederive the 
main results of the BCS theory (gap equation, thermody- 
namic potential, and collective modes), thus showing the 
ability of the 2PI RG scheme to access broken-symmetry 
phases. We find that the 2PI two-particle vertex ^^^^ is 
invariant under the RG transformation, while the fiow 
equation for the self-energy yields the BCS gap equa- 
tion. In Sec. IIVI we study one-dimensional (ID) sys- 
tems within the g-ology framework. At one-loop order, 
the 2PI scheme compares favorably to the RG scheme at 
high energy, but deteriorates at lower energy and even- 
tually breaks down. We identify the reason for this fail- 
ure and argue that the 2PI RG scheme can nevertheless 
be appHed to realistic quasi-lD systems like the organic 
conductors of the Bechgaard salt family. The conclusion 
is devoted to a discussion of the implementation of the 
2PI RG formalism beyond the examples considered in 
Sees, and ITVI 



RG approaches in the 2PI formalism have also been 
discussed in high-energy physics in a field-theoretical 
framework (see, for instance, Refs. I4^l4^48ll4fl( ) . 
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II. 2PI RG FORMALISM 

A. Luttinger-Ward functional and 2PI vertices 

We consider the partition function of a spin-i fermion 
system in the presence of external sources, 

Z[J] = j exp{-5[^] + ^VJi], (4) 

where the action S = Sq + 5int is defined by 



CKi ,a2 ,CK4 



i/'cr(r, t) if C = -, 
i/;;(r,r) if c = +. 



(6) 



Since the ^/''s anticommute, both Cap and Jq/3 can be 
chosen to be antisymmetric functions: Cap = —Cpa and 
Jap = — J pa- In the following, we denote by 



7 = {«,/?} 



(7) 



bosonic indices obtained from two fermionic indices a 
and (3. 

The single-particle Green function is given by the func- 
tional derivative of W[J] = lnZ[J], 



5W[J] 
(5J-V 



(8) 



Note that the definition of G differs by a minus sign from 
the usual definition in condensed-matter physics. The 
Legendre transform of W[J] is defined by 



T[G] = ~W[J] ~ iTr(JG), 



(9) 



where J = J[G] is obtained by inverting Eq. JSJ. To 
keep the notations simple, we shall denote J[G] by J in 
the following. Tr denotes the trace with respect to the 
fermionic indices, i.e. Tr(JG) = J2a p '^apGpa- r[G] 
satisfies the "equation of motion" 



6T[G] 

~sg7 



(10) 



as can be easily verified by a direct calculationi^ It is 
customary to write T[G] as 

r[G] = ^TrlnG- iTr(GG-i - 1) + 3'[G], (11) 

where the LW functional ^[G] is the sum of 2PI vacuum 
fiuctuation diagrams (or vacuum fiuctuation skeleton di- 
agrams), i.e. diagrams that cannot be separated into two 
disconnected pieces by cutting two lines. 

By differentiating Eq. ltTn)l with respect to the source 
J and using Eq. ^ fi^ we obtain 
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(12) 



with Ca0 the free propagator and Vaia2a3a4 the to- where 
tally antisymmetrized interaction vertex. ^ip'^Jip = 
^J2a p'^a'^apipp describes the coupling to the external 
bosonic sources Jap. The V'a's are Grassmann variables 
and the collective index a = (r, r, a, c) labels the posi- 
tion, imaginary time and spin projection along a given 
axis, as well as other possible internal degrees of free- 
dom. J2a = lo^T I '^'^^T.a.c "^here (3 = l/T is the 
inverse temperature and d the space dimension, c = ± is 
a charge index such that 



.(") 



^'^7l---7n 



5^"^T\G] 



SG. 



71 



■ 6G^ 



(13) 



are functional of G. Eq. lfT2ll defines a matrix multipli- 
cation with respect to the bosonic indices with / the unit 
matrix. Further relations between {W^^")} and {F^")} 
can be obtained by taking higher-order derivatives of the 
equation of motion ljl()|l . 

The 2PI vertices are defined by 



(n) 



(5(")$[G] 



(14) 



To order T^™, <i>(") is represented by all 2PI diagrams 
with n external (bosonic) legs 7i and 2m — n internal 
lines. These diagrams cannot be separated into two dis- 
connected pieces by cutting two internal lines (consider- 
ing every external leg 7,; — {a^, Pi} as a connected piece). 
The 2PI vertices satisfy the symmetry properties 



(n) 

7l---{ai,/3i}---7ii 



in) 

7l---{/3i,ai}---7„' 



in) 



(n) 



'7l---7i---7j---7„ — ■^7i---73---7i---7„- (15) 

The equation of motion lfT(1|l can be rewritten as a Dyson 
equation 



G^ G-^ J 'J 



■J71 



(16) 



with S-y — the self-energy. Note that the diagram- 
matic interpretation of $[G] as the sum of 2PI vacuum 
fiuctuation diagrams follows from Eq. II16|I (see, for in- 
stance, Ref. 1^3). Similarly, the equation v'^'^^W'^'^^ = I 
can be recast in the form 

^7172 ^ GaiP2Gpia2-G 

+ \ E G„,„3Gft&*?3^74<l2- (17) 
73,74 
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Eq. II17II is a Bethe-Salpeter equation relating the two- 
particle Green function W^"^^ to the 2PI vertex 

Equilibrium quantities are obtained for vanishing ex- 
ternal sources {J — 0). The equilibrium Green function 
G is determined by the stationary condition 



ST[G] 



5G^ 



0. 



(18) 



G 



It can be expressed in terms of the self-energy S = S|q 

$(1), 



(19) 



The thermodynamic potential is given by il 
-/3-MnZ[J = 0]=/3-ir[G]. 



B. 2PI RG equations 

We now take the free propagator C to depend on an 
infrared cutoff A which suppresses the low-energy degrees 
of freedom (|^k| < A), 

Ok 



C(k, zcj„) = — 



(20) 



C(k, iun) is the Fourier transform of C(r, r, — ; r', r', +), 
which we assume to be spin-rotation invariant, — 
Ek — M is the dispersion of the free fermions (with /i the 
chemical potential), and ujn a fermionic Matsubara fre- 
quency. 6k is a cutoff function such that 9k|A=o = 1 
and Ok|A=Ao = 0) where Aq = maxk |^k|- Physical quan- 
tities are obtained for A = 0, when all degrees of freedom 
are included in the partition function. W[J] and its Leg- 
endre transform T[G] now depend on A and satisfy flow 
equations as the cutoff A is varied between Aq and 0. 
Prom the definition of T[G] [Eq. lO], we deduce 



dl 



V[G] 



d_ 

dl 



W[J] 



dJ^ 5W[J] 
dl 5J^ 



1. 



-Tr 



dJ 
'dl 



G 



where I can be either A or a function of A, e.g. I — 
ln(Ao/A). Recall that J depends on G via Eq. ^ and 
is therefore a function of I. dW[J]/dl denotes the varia- 
tion of W [J] —InZ [J] which follows from the explicit A 
dependence of the propagator C, 



d_ 
dl 



W[J] 



1 



Z[J] 

--E 

Q,/3 



1 



-Tr(C-iG), 



(22) 



where the dot denotes a derivation with respect to I. 
imply 



dl 



T[G] 



1 



TriC-^G), 



(23) 



and, using Eq. Illl|l . 



(24) 



The LW functional is invariant under the RG transforma- 
tion. It is therefore a "universal" functional independent 
of the free propagator Ci^ This property has a simple 
diagrammatic interpretation. Being the sum of the 2PI 
graphs (with the internal lines corresponding to the vari- 
able G), <i>[G] depends on the interaction vertex V, but 
not on the non-interacting (A-dependent) propagator G. 

The thermodynamic potential satisfies the RG equa- 
tion 



1 d 
pJl 



1 

2P 



E 
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ST[G] 



5G~, 



G^ 



1 



Tr(G-iG), 



(25) 



where the last line is obtained using equations Ij23|l and 

Since <i>[G] is invariant under the RG transformation, 
the fiow of the_2PI vertices = ^^"^Ig is entirely due 
to the fiow of G, 



(n) 



d 
Jl 



(5(")$[G] 
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■SG, 



5("+i)$[G] 
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G 



^G 



(26) 



We thus obtain an infinite hierarchy of fiow equations 
for the 2PI vertices (Fig. HI . In all Feynman diagrams 
shown in this paper, a pair of neighboring external legs 
corresponds to a bosonic index 7i = {ai,Pi}. We shall 
always represent the 2PI two-particle vertex $7172 with 
(21) the two external legs 71 — {ai,(3i} (72 — {a2,/32}) on 



the left (right) hand side of the vertex: $7172 is 2PI as 
seen from left to right. 

Eqs. (I25l2(i|l should be supplemented with the initial 
values of the thermodynamic potential and the 2PI ver- 
tices at A = Aq. Since the cutoff function 9k is chosen 
such that C\ao = 0, we easily deduce from the diagram- 
matic expansion of the LW functional that S|ao = 
(i.e. G|ao = 0), ¥^^\ao = V, and $(")|ao = for 
n > 3. Using Eq. lfTT|l . we also obtain fijAo = ^^oIaq 
where flo = (2/3)~^TrlnG is the non-interacting thermo- 
dynamic potential obtained from the action Sq. To avoid 
computing flo, we shall always calculate AD, = V, — Qq. 

For practical calculations, one has to truncate the hier- 
archy of fiow equations Ij26|l by retaining a finite number 
of low-order vertices. The simplest non-trivial truncation 
is discussed in the next section. Unless mentioned other- 
wise, we now drop the "bar" above Green functions and 
vertices since we shall only consider the case J = 0. 




FIG. 1: Diagrammatic representation of the hierarchy l|26|l of 
flow equations satisfied by the 2PI vertices "1> 



(n) 



Deriva- 



tion with respect to I are indicated by dots (vertices) and 
slashed lines (propagators). A pair of neighboring external 
legs corresponds to 7i = {at, Pi}. 








FIG. 2: Diagrammatic representation of as a function 
of within the one-loop approximation. shown as a 

black dot, is 2PI as seen from left to right. Signs and symme- 
try factors are not indicated. 



C. One-loop RG equations 

1. Flow equations for the 2PI vertices 

One-loop RG equations are obtained by neglecting 
for n> A. This reduces the hierarchy of equations 
(Eel to 



t = iv$(^) G 

^71 2 7172^72 

72 

$(2) ^ iv$(3' G 

7172 9 / J 717273 



In order to close this system of equations, we need an 
approximate expression of '^^^^ in terms of and S. 
Let us start with the second-order contribution to the 
LW functional, 

^[G] = — — ^ VQi;33Q4/32Va3/3ia2/34G'7l ^72^*73 G^4- 

71,7273,74 



By taking the third-order functional derivative with re- 
spect to G, we obtain 



oneself that the overall factor —1/2 in Eq. (I29II remains 
unchanged when V is replaced by ^'•^■'Iiioop is 

shown diagrammatically in Fig. [2l 

Diagrams contributing to <i>[G], S, and up to 
third order in the bare interaction amplitude V are shown 
in Fig. El The 0{V'^) contribution to <i>(^) is included in 
the one-loop approximation, but among the three dia- 
grammatic contributions to order only the first one 
is retained. The other two are not of the form the 
second one involves a two-particle-reducible two-particle 
vertex, while the third one is clearly not of the required 
type. Thus, a given diagram contributing to the LW 
functional will generate diagrams for which may or 
may not be included in the one-loop approximation. The 
latter is therefore not a ^-derivable approximation as it 
is not based on a truncation of the LW functional. A 
detailed study of conservation laws and Ward identities 
in the 2PI RG scheme is beyond the scope of this work 
and remains to be done. 

It should also be noticed that the one-loop approxi- 
mation, like any approximation of the 2PI vertex 
leads to a violation of the crossing symmetries of the 
(28) two-particle Green function M/^^) (g_g^ W^if^j.^./j^ 7^ 



(27) 



(3) 

717273 



Q ^ ^ G-y4 [Vq,iQ!3Q!2Q4 ^1 



/33/32/34 



^lQ3/32a4 ^1/33112/34 ^1/34/3203 y(3ict4a2/33 
-V"ai/34a2a3V"/3iQ4/32/33 " (^S ^ Pi)]- (29) 

Replacing V by in Eq. (I29L ^^ we obtain the one-loop 
approximation of 



$(3) 
7172 73 



Hoop 2^ ^^L"^"^"^" 

74 



(2) 



2a4^/3i/33/32/34 



(2) 



(2) 



0103/32114 /3i/33a2/34 



(2) 



(2) 



(2) 



(2) 



ai/34a2a3 /3ia4/3; 



ai/34/32Q3 I3iona203 
2/33 - ("3 ^ /33)] • 



(30) 



By working out the symmetry factors of various diagram- 
matic contributions to C^*^^^ and <i>'^'^\ one can convince 



-W, 



(2) 



) and the IPI two-particle vertex. 



From Eqs. 112713011 . we finally obtain the one-loop RG 
equations 



t - iV$(2) G 

71 — n / .^7172 '-^72; 



$(2) 
7172 



2 

72 

2 ^ 

73,74 



(G73 G^^ 



G73 G~f^ ) 



L Qi 



(2) 



0304/32 '^/3i/33/34a2 



(2) 



(2) 



aiQ3a4Q2 /3i/33/34/3' 



(31) 



Eqs. l(3T|l are shown diagrammatically in Fig. ^ There 
are two differences with respect to the one-loop RG equa- 
tions obtained within the IPI RG schemetS2iS4 (i) the flow 
equation for the self-energy involves the 2PI two-particle 
vertex and the derivative G of the Green function in- 
stead of the IPI vertex and the "single-scale propaga- 
tor" S = —GC~^G; (ii) the one-loop contribution which 
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FIG. 3: Diagrams contributing to $[G], E-yj, $V72 

(3) 

72 73 up to third order in the bare interaction vertex V 
(shown as a (small) black dot). Diagrams obtained by ex- 
changing external legs {oi ^ are not shown. 




FIG. 4: One-loop flow equations for and $7172. Here 
and in the following figures, diagrams obtained by exchanging 
the slashed and non-slashed lines in the one-loop diagrams 
contributing to "t'^' are not shown. 

would give a two-particle- reducible contribution to 
is absent. 

The procedure we have followed to obtain the one- 
loop RG equation is not unique. Owing to the anti- 
symmetry of V, one could write -Vai/33a4/32^a3/34a2/3i 

or 144/33ai/32^Q3/34Q2/3i instead of 14i^3a4/32Ka;3/3ia2/34 in 
Eq. Ij28|l . This explains while the one- loop RG equation 
obtained in Ref. |lj 

$(2) _ V G G f$^^) 

7172 ~ ^73^74 L^a4a3ai/32 /34/33/3ia2 

73,74 

~ ^a4a3aia2*^k/33/3i/32] ' (^^2) 

differs from ours. This equation can be represented di- 
agrammatic as in Fig. 0) but with the vertices in 
the one-loop corrections being 2PI as seen from top to 
bottom. Eqs. II31I32|I lead to different diagram resum- 
mations. While the two 0(V^) contributions to in 
Fig. El are generated by the equations iPIJ; the first one 



is not if one uses Eq. More generally, similar dia- 

grams with an arbitrary number of loops are not included 
in Eq. ll32|l . These diagrams play a crucial role in most 
applications of the RG approach to interacting fermion 
systems. For instance, they describe the exchange of spin 
fluctuations in a conductor with short-range antiferro- 
magnetic order and may lead to d-wave or other types of 
unconventional superconductivity. For this reason, we do 
not expect the one-loop approximation based on Eq. \,V1\ 
to give reliable results. 

2. Relation to the IPI RG scheme 

In this section we show how, starting from Eqs. (EH, 
we can reproduce the one-loop RG equations for the 
IPI two-particle vertex 7'-'*-' obtained within the IPI RG 
scheme. 7^**^ is deflned by 

^'fnil2 = ~ 51 ^"103 G'^l/33Ga2a4G'/32/34 77374; (33) 
73,74 

where 

'^c^7l72 = {'^a^^^3li^a2i^h) - GatPiGa^p^ 

— ^7172 ~ G 01113^0(3^012 + Gotioi2Gl3iP2 (34) 

is the (fully) connected two-particle Green function. The 
Bethe-Salpeter equation Ijl7|l implies that G^f^ satisfies 

2 [^7174 ^ 2 ^ ^^7173 '''7374] >-^c,7472 ~ 

74 73 

~ ^ 5^ -^71 73*^*73 74 1^74 72; (35) 
73,74 

where 

117172 = Gail32Gf3ia2 ~ GqjQ2 '-'/Sift • (36) 

In matrix form, Eqs. H33I35|I read 

G(4) = -n7(4)n, 
(/ + n$(2))G'(4) = -n$(2)n. (37) 

We therefore obtain 

7(4) ^ ci>(2)(/ + na>(2))-\ 

$(2) = 7(4)(/_n7W)"\ (38) 

From these two equations, we deduce the following fiow 
equation for the IPI vertex 7^"'-', 

^(4) = _^(4)n^(4) _^ ^(4) ($(2)^-1^(2) ^^(2)^-1^(4) 

= -^(4)n^(4) + (/-^(4)n)$(2)(j_n7(4)).(39) 

We can now reproduce the one-loop RG equations de- 
rived within the IPI RG scheme by expanding the rhs 



7 



of Eq. to second order in 7^"'^ Since <I>(^) = 

7^^) + 0[(7^^^)'] and $(2) = 0[($(2))2]^ gq. ^ gives 

^(4) ^ _^(4)n^(4) _^ ^(2) ^j^g^g^ Q^^g^ (^(4))2^ ^(2) jg 

given by Eq. with <i>(2) replaced by 7'^''-' in the rhs. 
This eventually gives 



'7172 



2 $Z (^"3/34^/3304 + Gag/j^G/jjQ^) 

73,74 

L 'Qi/3ia3/33 'a4/94a2/32 'aia^psp^ ' pi [3 4,014,01-2 

T/3i/34a4/32J ■ 

(40) 



To lowest order in 7^^\ the flow equation (|3Tll for the 
self-energy becomes 



(41) 



Eqs. Ij4()l41|l agrees with the equations derived within the 
IPI RG schem922*2i with the "single-scale propagator" 
S = -GC-^G replaced by G. 



D. One-loop RG equations in the normal phase 

In this section, we consider the one-loop RG equa- 
tions in the absence of broken symmetry. We de- 
note the position, time and spin indices by X so that 
a = (X,c), = ^*{X), ^X^) = V(^), and 

^ dX — Jq dr J d^rJ2a- ^he single-particle Green func- 
tion is then given by 

G{Xi-,X2+) = mXi)r{X2)) = G{Xi,X2), (42) 

with G{Xi+,X2-) = -G{X2,Xi). G{Xic,X2c) van- 
ishes in the normal phase. G satisfles the Dyson equation 

GiXuX2) = C{Xi,X2) 

-J dXzdXiC{Xi,X:i)^{X:i,Xi)G{Xi,X2) (43) 

with 

E(Xi,X2) = S(Xi+,X2-) (44) 



the self-energy. In the normal phase, the 2PI two-particle 
vertex ^'^^^ (XiCi, X2C2, X3C3, X4C4) vanishes if ^ 
0. Distinguishing between the particle-particle (pp) and 
particle-hole (ph) channels, we introduce 

$(,2p)(Xi,X2,X3,X4) = '^^^\Xi+,X2+,X:,-,Xi~), 

3>;,'J(Xi,X2,X3,X4) = <^'^-'\X^+,X2-,X^ + ,X^-). 

(45) 



Note that in Eqs. l|l5)l . we have singled out one of the 

two ph channels. The 2PI vertex in the other ph channel 

(2) 

is related to <i>pj^ by 



$^"^(Xi-f,X2-,X3-,X4 + ) = -$|,'i](Xi,X2,X4,X3). 

(46) 

Since the one-loop approximation conserves the cross- 



ing symmetry ^^l^^^p^} = "^Tiifeas} (^^ obvious from 
Fig- El) it is possible to consider a single ph channel. The 
symmetry properties of $7^72 [Eqs. Ill5|l ] imply 

<I>(2)(Xi,X2,X3,X4) = -$(2)(X2,Xi,X3,X4) 



(2) 



$;,'j^)(Xi,X2,X3,X4) 



-<i>(2p)(Xi,X2,X4,X3), 
$2,^(^3, X4,Xi,X2). (47) 



The RG equation for the self-energy can be written as 



t{Xi,X2) = 



f dXsdX^GiXscXic) 

X$(2)(J^^ + ^X2-,X3C,X4C) 

- J dXsdXiGiXi^Xs) 

X<i>('h^(Xi,X2,X3,X4), 



(48) 



where c — ~c. A similar calculation for the two-particle 

(2) (2) 
vertices <i>pp and ^^^^ yields 



<i>(2p)(XO = -i I dX[dX'2dX:,dX',[GiX!2,X[)GiX',,X',) + iG^G)][^'^^{X^,X[,X',,X^)^^^^ 

+$g(Xi,x^,x^,X4)<i>(2p)(X2,x;,x^,X3) - (X3 ^ X4)], 

- -I JdX[dX:,dX!,dX',[G{X',,X[)G{X',,X',) + iG^G)][^^^^{X,,X[,X!,,X 

+^'^^{X,,X!„X'„ X,)^^^{X[,X2,X,, X'^ + $g,'(Xi, X^, X3, X^)$(1'(X( , X2, XI X,)] , (49) 



(a) 



(b) 



> ■ X) 





•IT- 




(0 





FIG. 5: One-loop RG equations for E (a), $pp (b) and ^ J (c) in the normal phase. 



where we use the short-hand notation = $(2) (^^^^2, X3, X4). Eqs. ^ are shown diagrammatically in 
Fig.H 

In spin-rotation invariant systems, it is convenient to write <i>pp and in the triplet/singlet and charge/spin 
basis, respectively. 



(72 .(74 ■) 



(50) 



{X = (XjCr), X — (r,T)), where 



(51) 



and cr = (a^, a^, a^) stands for the Pauli matrices. From Eqs. ifTTjl . we deduce the following symmetry properties 



^f^{x-Y,Xi,X'i,X^) 

*ch!sp(2;i'2;2,a;3,a;4) 



-$J: '^(x2,xi,2;3,a;4) = ^(xi, 2:2, X4, 2:3), 

^f\x2,Xx,Xz,X^) = ^f^{xx,X2,X^,X-i)-, 



*ch,sp(2;3,a;4,a;i,a;2). 



(52) 



Performing the sum over spin indices in Eqs. II48I49|I . we obtain the flow equations satisfled by S, and ^^h^gp. 



S(a;i,a;2) = -2 / dx'idx/^G(xt,,X'i)^'^^^(xx,X2,xz.,x^), 



(53) 
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where 



<i>f ( 



= — - J dx\dx'2dx'^dx'^B{x'2,x\,x'^^x\) 

X ) {x,,x'„x',, x,)<^^2 4, ^3) + 2$f )$g) + ci,(2)$g) 

~ J dx -ydx '^dx "^dx ^ Bi^X'^ ^ x-^^ "^37 •^4. ) 
X x'l, 4, a:4)<i>fh^(x2, 4, 2;^, ^^s) + 3$|'^$g)) 

- ~ ^ dx -^dx ^dx ^dx ^ B (^x^ , 7 "-^3 : •^4 ) 

X (^^^[^\xux[,x'3,xM^\x'^,X3,X2, X'2) + ^<i>(')$f 

+ ($ (xi , a;'2 , , X4 ) $fj (x'l , X2 , X3 , x^, ) + 3ci>g) <i>g) 
^ip (a^i) ^ ^\ I dx[dx2dx3dx'4B{x2,x[,x'3,x'^) 



i$g) (xi, x'l, 4, a;4)<i>r (4, =^3, 0:2, 4) + \<^?^^P + l^P^f^ 
+ (*ch (2^1 ' 2;'2 , 2:4 , X4 ) $g) (a; 1 , X2 , X3 , ) + 

*^sp^*^ch ^" 2$g''$g^ 

+ [^^^{x^,x'2,x,, 4)$g) (x'l, X2, 4, ^4) + 1>g)<h^ - 2$g)$g) 
i3(a::i,X2,2;3,a;4) = G'(a;i, X2)G'(a;3, X4) + (G ^ G). 



(54) 



(55) 



Due to spin-rotation invariance, the single-particle Green function Gaia2{xi,X2) — S„^^„2G{xi,X2) and the self-energy 
Scr^o-2(3;i,a;2) = Sai,cr2^{xi,X2). In Eqs. Ij54|l . we have grouped inside parenthesis terms with identical dependences 
on the variables Xi,x^ and thus avoided unnecessary repetitions of the latter. 



E. Response functions in the normal phase 

Once we have solved the RG equations H53I54|I and determined the self-energy E and the 2PI vertex we can 
obtain the two-particle Green functions VF^^^ and the response functions from the Bethe-Salpeter equation l|17|l . 
We define triplet- and singlet-pairing fields, and charge- and spin-density fields by 

{'^t(xi)iIjt{x2) if 1^ = 1, 

^EaMxi)Mx2) if 1^ = 0, 
'>Plixi)'>Pl{x2) V = -1, 



--i=^cr'0^(xi)V'ff(a;2), 

p(xi,2;2) = ^^l{xi)^a{x2), 

S(X1,2:2) = ^ ll)l^{xi)cT^^^„^ll)„^{x2), (56) 
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(a = — cr), and the corresponding response functions 

W^pp!iiii(^0 if '^--l, 



Cl ,(73 

Xsp(a;z) = (5''(a;i,a;2)5''(x3,a;4)) = ^ ^^i,<T2'^J^3,'T4^ph!cT2<Ti^4CT3 (2^2, 0:4, X3), (57) 



CJ"! :Cr2 if^S iCr4 



where H/j^pV,(a;«) = M^^, Ha^i", 2:2-, 2:3+, a;4+) and W^^^^^^^ (x,) = M^^,^(a;i- 

, a;2+, xy,—,Xi+) are the two-particle Green 
functions in the pp and ph channels, respectively. They satisfy the Bethe-Salpeter equations [see Eq. ifTTjl ] 

Appier; i-Xi) ^(T 1 ,CT4 ^fT2 ,Cr3 

G(a;i,a;4)G'(a;2,a;3) - 5 G(xi,a::3)G'(x2,a:4) 
+ ^2;W4^^4'^a;4G(xi,a;i)G(a;2,a;2)^pp,^i<r2aa' (a^i, 4. 4. 2:4)W^pp!<T<T'<T3a4 (4. 2^4. a;3, 2:4), 



+ dx[dx',dx',dx\G{x,,x[)G{x',,X2)^^^l,^,^,,, ix[,x'2,x',, 4)W^S.v.3-4 (4, 2^:^x3, X4). (58) 

I7,(t' 

From Eqs. (I58II we deduce 

Xti=t,six.i) = G{xi,X4)G{x2,X3)tG{xi,X3)G{x2,X4:) 

J dx[dx2dx'2dx'4^G{xi, x[)G{x2, X2)^^^\x[, X2, X3, x'4)Xt,{x4, X3, X3, X4), 
X/i=ch,sp(a;j) = -2G(a;4,a;i)G(a;2,a;3) 

+2 / dx[dx2dx3dx'4^G{x2,x[)G{x2,xi)^^^\x[,X2,X3,x'^)xp,(x';j,X4,X3,X4)- (59) 



Eqs. lf59jl enable ones to determine the response functions and considering this equation in the pp and ph channels, 
from the knowledge of the single-particle Green function we obtain 
G and the 2PI vertex $(2) . 



ll=t,si'^^) = '^^i^\x^)-^ [ dx'^dx'2dx'3dx'^ 



2 . 

x$|f)(a;i, a;2, 4, 2;'2)G(x'i, 4)G(x^, a;;) 

In many cases, useful information can also be drawn y.-f^if^ {x'^,x''^,X3,X4), 

from the IPI two-particle vertex 7^"^^ Rewriting the re- r 
lation ^ between 7''^) and $(2) %=ch,sp{x^) = '^^i^'i^i) + 2 J dx[dx2dx3dx'4 

x$|f)(a;i, a;2, 4, 4)G(xj5, a;'i)G(x^, s^) 
X7u'*H2;4,4,2;3,a;4). (61) 



'^7172 ~ ^7172 2 ^7i73^a3/34G'/33a4 77472' 

7374 



F. Ginzburg-Landau expansion 

An essential feature of the 2PI scheme is the possibil- 
ity to continue the RG flow in a broken-symmetry phase. 
This will be illustrated in the next section in the frame- 
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work of the BCS theory, and further discussed in the 
concluding section. As a byproduct, one can also derive 
the Ginzburg-Landau expansion of the thermodynamical 
potential in the vicinity of a phase transition. The inter- 
est of such an approach is that only the solution of the 
RG equations with no symmetry breaking is necessary. 
When several instabilities compete at low temperature, 
the derivation of the Ginzburg-Landau expansion is ex- 
pected to be much simpler than the full solution of the 
RG equations in the broken-symmetry phase. 

We introduce the Legendre transform F[T,] of the LW 
functional ^[G], 



In the vicinity of the phase transition, where A is small, 
we expand r[A] to fourth order. 



r[A] = Tn 
1 

~6 



-— V A A 

2! 22 ^ -^7172^71^72 
7li72 



Tr(G^A)3 



1 1 
3! 23 



'^7?72 73^7l^72^73 



-hiTV(G^.Ar + l^ 



71:72,73:74 



P(4) 
7172 7374 



X A-j-j A^2 A^3 A^^ , 



Fp] - a>[G] + iTr(GS), 



(62) 



(67) 

where I3~^Tn = P~^T[E,n] = ^In is the thermody- 
namic potential and Gat the Green function in the normal 
phase, and 



where G = G[E] is a functional of S obtained by inverting 
= 5^[G]/SG^. Here and in the following, we assume 
vanishing external sources {J = 0), and denote G, E by 
G, E, etc. F[E] satisfies 



in) 



(5(")F[E] 



(5E^j • • • (5E^^ 



(68) 



sm 

(5E-V 



-G, 



(63) 



The functional F[E] allows us to rewrite r[G] as a func- 
tional of the self-energy;^ 



If, as in most one-loop approximations, one ignores the 
normal phase self-energy (Gat = G), the evaluation of 
the terms Tr(GArA)" does not raise any difficulty. The 
determination of r[A] then requires the calculation of 
the coefficients i^'"). These can be related to the 2PI 
vertices by taking functional derivatives of Eq. Ij63|l 
with respect to G. The first-order derivative gives 



(5(2)$[G] 



r[E] 



iTrln(G-i + E) + iTr(GE) + $[G] 
1 



2 ^3 SGj^SG^^ 



(5(2)i^[E] 



7172 J 



-Trln(G-i+E) + F[E], 



(64) 



I.e. 



\ / 71-; 



7172 



.$(2) 
7172 ■ 



(69) 



(70) 



which is stationary at the equilibrium self-energy E = E, Higher-order derivatives yield 



6m 



= 0. 



(65) 



7172 73 



P(4) 
71727374 



- y a>(^ 

S /-^ 7i 



(3) ^(2) ^(2) p{2) 
7172 73 7171 72 72 73 73' 



7i:72:73 



We now write the self-energy E = Ejv + A as the sum 
of a normal part Eat and an anomalous part A which 
violates some symmetries of the normal phase. A is an 
order parameter for the phase transition. It can include 
different kinds of symmetry breaking, such as antiferro- 
magnetism and superconductivity. If we are able to solve 
the flow equations in the normal phase, i.e. flxing A = 0, 
we can consider r[A] = r[Ejv -t- A] as a functional of the 
(unknown) anomalous self-energy A. The latter is then 
determined from the stationary condition H65II . 



1 

16 



^(4) ^(2) ^(2) ^(2) ^(2) 

7172 7374 7l7l 72 72 7373 74 74 



E 

7;, 72:73.71 
S 7172 73 L 747171 72 72 73 



(2) 

73 



7i:72.73 



5r[A] 

5A^ 



= 0. 



(66) 



A crucial point here is that the 2PI flow equations with 
A = can be continued below the actual phase transition 
temperature Tc, since the 2PI vertices do not become sin- 
gular at the transition. The functional r[A] can therefore 
be determined below T^, where the stationary value of A 
[Eq. (jSEl] is flnite. 



,^(2) ^(3) ^(2) ^(2) p{2) pi3) 1 
I'lll 7472 72 7373 7171 72 72 7473 73^ ' 

(71) 

These equations simplify within the one-loop approxi- 
mation, since <i>*^"^ vanishes for n > 4 and ^'-^^ can be 
expressed in terms of F'^^) a^jj ^(4) g^j.g ^j^g^ essen- 
tially determined by ^'^^^ and F*^^' . Writing the equation 
(j)(2)^(2) _ _/ ^jjg pp g^jjj pjj channels, we obtain 

i I dXdX'<^>^^^{XuX2,X,X')Fl;l\X',X,X3,X,) = 
- 6iXi - Xi)S{X2 - X3) + S{Xi - X3)S{X2 - Xi), 
J dXdX'^^^^l iXi,X2,X, X')F^^ {X', X, X3, Xi) = 

- 5{Xi - Xi)5{X2 - X3), (72) 
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where 



FllHX,,X2,X3,Xi) = F(2)(Xi-,X2-,X3 + ,X4 + ), 
F^^{Xi,X2,X3,Xi) = F(2)(Xi-,X2 + ,X3-,X4 + ). 

(73) 

For a spin-rotation invariant system, we finally deduce 

dxdx'^l^^ {xi,X2, X, x')fI:'^^ {x' jX^X^^Xi) = 

2[5{xi - Xi)5{x2 - X3) - 5{xi - X3,)5{x2 - x^)], 

dxdx'<^^^^ {xi,X2,x, x')F^'^\x' , a;, 0:3, a;4) = 
2[5{xi - xa)5{x2 - 2:3) + ~ x^)6{x2 - x^)], 
dxdx'^''^ {xi,X2,x, x')F^y^ {x' , x,xy,,Xi) = 



- ■^5{xi ~ X4)5{X2 - X3), 
dxdx'^fj (Xi, X2, X, x')i^ip^ {x\ X, X'i, Xi) = 

1 




FIG. 6: Luttinger-Ward functional $bcs[G] in the BCS the- 
ory. 



Alternatively, one can start directly from the LW func- 
tional (Fig. El 

$Bcs[G]==A J dxGTi{x+,x+)Gn{x-,x-), (77) 

where 

Gn{x+,x+) = {■4;*{x)t(;l{x')), 
Gn{x-,x'-) = {^lix)^P^{x')). (78) 

One then obtain ^^1^273 = 0; which leads to Eq. Ij76|l . 
In the case of a local interaction, the only non- 

/2) 

vanishing part of $pp reads 



Six, ~ x,)S{x2 ~ X,). (74) $(2) (^J ^ m.Mx, ~ X2)5ix2 - x,)Six, - X,) 



(2) (2) 

where F^ ,,' are the triplet and singlet parts of i^pp , and 

-^ch sp charge and spin parts of • In the next 
section, we shall use Eq. (|67ll to reproduce the Ginzburg- 
Landau expansion of the thermodynamic potential of a 
BCS superconductor. 



III. BCS THEORY 

The aim of this section is to reproduce the main re- 
sults of the BCS theory using the 2PI RG equations. We 
consider a 3D system described by the action 



^pp,^;^^^ -pp 

_ _^(2) _ _i(2) _ ? (2) _ . 

(79) 

i.e. I>s^^ = 2A and "if ^ = 0. The self-energy has two 
non-vanishing elements, S||(a:+,j/+) = — I]j^|(?/+, x-f ) 
and Tti^{x—,y~~) = — S|| (y— , ), determined by the 
RG equations 

^ni^+^y+) = ^^{x ~ y)Gii:{x-,y-), 
t^^ix-,y-) = XS{x-y)Gtiix+,y+). (80) 

Defining the superconducting order parameter A by 



S 



dx ^i)l{x) ( 9^ - /i 



2m 



+A / dx'il:*^{x)ipl{x)i:i{x)i:^{x), 



(75) 



where m is the fermion mass, /i the chemical potential, 
and A < the amplitude of the local attractive interac- 
tion. This singular interaction is regularized by means 
of an ultraviolet cutoff acting on the fermion dispersion: 
ICkl < Ao, a = kV2m-^. 



S|j,(.T+,y+) = 5{x-y)A, 
SxT(-^-,y-) = 5{x-y)A\ 

the fiow equation becomes 

K^\Gi^{x~,x-). 



(81) 



(82) 



Integrating Eq. Ij82|l between Aq and A, we obtain the 
gap equation 



A 



(83) 



A. BCS gap equation 

There are two equivalent ways to derive the RG equa- 
tions for the BCS theory. One can start from the one-loop 
equations and neglect the ph channel. Since the one-loop 
contribution to ^^p* involves only the ph channel, it van- 2x2 matrix in reciprocal space. 



where V is the volume of the system, and k — (k, ia;„) 
with w„ a fermionic Matsubara frequency. 

We use the standard Nambu notation, "^{x) = 
(■!/'|(x),-!/'|(a;))"^, to write the inverse Green function as a 



ishes in the BCS approximation, i.e. 

$(2) = 0. 
^pp 



(76) 



C-\k) 



-C-\-k) J ' 



(84) 
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where C(k) is the free propagator defined by Eq. Il2()ll . 
For intermediate calculations, it is convenient to assume 
that &k is a smooth cutofi^ function which does not vanish 
(©k 7^ 0). Final results depend only on C(fc) (and not 
C~^{k)) and are well defined even for a hard cutoff Ok — 
0(ICk| — A) (6 is the step function). Inverting G'^^ , we 
obtain 



Gaa {k , k- 



GiT(fc+,fc+) = eK 



El' 



At 



where 



El' 



(85) 



(86) 



Using Eqs. (I85II . we can rewrite the gap equation Ij83|l as 



A 

V ■ 



2 



(87) 



When A 7^ 0, this equation becomes 



1 



= tanh 



lAI 



where we have taken 0k = 0(|Ck| ~ A). The density of 
states in the normal phase N{^) = yj^k^i^ ~ '?k) has 
been approximated by its value iV(0) at the Fermi level, 
and A — N{0)X is a dimensionless interaction constant. 
Eq. Ij88|l can be solved exactly at T = 0, 



|A(T = 0)| = |Ao| 1 



2A 



1/2 



e 



|Ao| 



A 



(89) 



2Aoei/^ is the 



assuming |A| < Aq. |Ao| = |A|a=o 
zero-temperature BCS gap. 

When A(T, A) = 0, the singlet response function 
Xs{x , X , y , y) defined in Sec. Ill El is given by the Bethe- 
Salpeter equation l|59|l . In Fourier space. 



Xsiq) 



2n-+(g) 

l + An-+(g) 



(90) 



{q — i'ij'i-^v) with uji, a bosonic Matsubara frequency), 
where 11 (g) is the Fourier transform of the pp propa- 
gator 



U-+{x-y) = [C{x,y)f 



(91) 



The condition for Xs{q — 0) to diverge, 1 + AH '^{q = 
0) = 0, is equivalent to the linearized gap equation 
[Eq. IHHl with A ^ 0+]. 

Fig. 13 shows the superconducting order parameter A 
and the singlet response function Xs{q = 0) (in the regime 
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FIG. 7: Order parameter A and singlet response function 
Xs (for A(r,A) = 0) versus A for T = (solid lines) and 
T = 0.9Tc (dashed lines). 



where A(T, A) = 0) at T = and T = O.QTc where is 
the superconducting transition temperature. Xs diverges 
at the threshold value Ac(r) of the cutoff below which 
the gap A becomes finite and reaches the BCS value for 
A = 0. The transition temperature is determined by 
Ac(T'c) = 0. As the 2PI RG equations determine A and 
they are not plagued with divergencies and can be 
continued down to A = for any temperature. Broken 
symmetry is signaled by a finite value of the anomalous 
self-energy A below Ac(T). The divergence of the singlet 
response function Xs is seen only when the 2PI vertex 
^^"^^ is fed into the Bethe-Salpeter equation relating Xs 

to 4p. 

Note that we have picked up the nonzero solution of 
(IsTIl by hand. Alternatively - and this is how one should 
proceed in more complicated situations - one can directly 
solve the fiow equation Il82ll i^ 



A 

X 



A 



VA2 



= tanh (^^VA2 + A2^ 



AA^ 



E^ 



^tanhf^i;)- 
E \2 J 2cosh' {(3E/2) 



(92) 



(E = ^y ^'^ + A'^ and we assume the gap A to be 
real), together with a symmetry-breaking initial condi- 
tion A(Ao) = e. For e/A(A = 0) < 10"^ the solution 
of H92II cannot be distinguished from that of H88|l when 
plotted on the same graph. A larger value of e leads to a 
smearing of the singularity at Ac(T). 



B. Thermodynamic potential 

We could determine the thermodynamic potential = 
l3^^r directly from Eq. JTIJ using the expression Ij77|l of 
'i'BCS- To illustrate how the 2PI RG scheme works, we 
shall instead use the flow equation (l25|l and the Ginzburg- 
Landau expansion ijfiTjl of the thermodynamic potential 
near a phase transition. 
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1. Ground state condensation energy 

The ground state condensation energy AQ = fl — 
satisfies the RG equation [see Eq. lf25| ] 



An = -—tt[c-\g-c)] 

2 ^ lApekOk 



2 

n 



(93) 



where the last line is obtained for T = 0. A direct evalu- 
ation of this equation by replacing the sum over k by an 
integral over ^k leads to ambiguities because of factors 
6(0). We therefore rewrite Eq. l|93jl as 



A(i = -dAj2Ja + \A\^e 



)2 



(94) 



where 9a denotes a derivation with respect to A at fixed 
|A|. The calculation is then straightforward and yields 



An 



~2VN{0) A - VA2 + |Aj2 
-2F7V(0)(2A- |Ao|)e' 



A , (95) 



where we have used Eq. l(89|l . Integrating this equation 
between Aq and A, we finally obtain 



An = -2VN{0)(^-^ 



A 6 



|Ao| 



A 



(96) 



For A = 0, we recover the expression of the condensation 
energy An ~ ~^VN{0)\Al\ in a BCS superconductor. 



2. Ginzburg-Landau expansion 

In the vicinity of the superconducting transition, the 
thermodynamic potential can be calculated using the 
general method discussed in Sec. Ill Fl This approach re- 
lies on the 2PI RG equations in the normal phase (A = 0) 
continued below the actual transition temperature Tc (i.e. 
with S = Eat = and ^^^^ = const in the BCS model). 
The anomalous self-energy A^^^^ is defined by 

Acria2ixiCl,X2C2) = S{xi - X2)6cri ^ai^ci ,C2 

xai{S,,,+A-5c,,-A*). (97) 



Using H97|l . one finds 

Tr(GwA)2 = 

Tr(GwA)4 = 

- V a a ' 



-4|Ap^C(fc)C(-fc) 

k 

-41/^1 ApiV(O) In 
4|A|4^[C(fc)C(-fc)]2 



VP\A 
lAP 



,4 7C(3)A^(0) 



-Vf3- 



lAI 



(98) 



A 



where 7 ~ 1.78 is the exponential of the Euler constant 
and C{z) the Riemann zeta function (C(3) ~ 1.2). To 
obtain the last fine, we have used the relation iTtHi be- 
tween and fP and Eq. From Eqs. II67I98|I . we 
recover the Ginzburg-Landau expansion in a BCS super- 
conductor, 



An 



|A| 



iV(0)ln 



27A0 
ttT 



7C(3)iV(0) 4 
167r2r2 I I ■ 

(99) 



C. Response function and collective modes in the 
superconducting phase 

For A < Ac(T), where A(T,A) ^ 0, the singlet super- 
conducting response function is defined by 



Xs 



2 



a\a2 



CTl ,0-2 

(2) 



(Ti<Ti<T2<T2 (2^1 ' ^1 ' ^1 ' ci , a;2, C2, a;2, C2). 

(100) 

Using the Bethe-Salpeter equation ifTTjl . we obtain 

X-+(g) - 2n-+(q)-An"+(q)x-+(g) 

-An— (g)x^+(g), 
x++(g) = 2n++(g)-An+-(g)x++((z) 

-An++(g)x7+(«), (101) 

where W'^'^'^^q) is the Fourier transform of 

n-+(x-y) = n+-(j/-x) =GTT(x,2;)Gii(x,2/), 
n (x-y) = -G4t(x-,y-)Gtx(a;-,?;-), 
n++(x-y) = -GTx(.x+,y+)GxT(a;+,2;+). (102) 

Collective modes are obtained from the poles of X^^^^('7)) 

[1 + An-+(q)][i + An+-(q)] - A2n++(g)n~(q) = o. 

(103) 
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Using 

n-+(g = 0) = ll+-{q = 0) 
U—{q = 0) = n++(q==0)* 

one verifies that Eq. 1110311 is satisfied for q = and 
Wi, = 0. For any finite value of the gap A (i.e. for 
A < Ac), we therefore obtain a gapless (Goldstone) mode 
(Anderson-BogoUubov mod o^^'^^ ). Note that the Gold- 
stone theorem is ensured by the fact that the BCS theory 
is a ^-derivable approximation. 




FIG. 8: One-loop RG equations for the 2PI vertex "3?'^ in 
the Cooper (a), Peierls (b) and Landau (c) channels in a ID 
system. Solid and dashed lines indicate right- and left-moving 
fermions, respectively. [In the pp loops appearing in (b) and 
(c), a sum over the left/right indices is implied.] 



IV. ONE-DIMENSIONAL SYSTEMS 

In this section, we consider a ID system with the action 
(for a review on ID systems, see Ref. 12^ 

k.r,a 

Sint = ^ {92Sr,r' + giSr,f') 

xrrAk + qWf.'{k' - g)V'rV'(fc')VVa(fc|J05) 

where L is the length of the system, k — (fc||,iw„), 
and q = (giijiwiy). k\\ and denote momenta, w„ and 
ujy fermionic and bosonic Matsubara frequencies, respec- 
tively. The index r distinguishes between right (r = +) 
and left (r = — ) moving fermions. Cr(fc||) — ^r{k\\) — M = 
VF{Tk\\ — kp) is the dispersion law, Hnearized around the 
two Fermi points ±/c_f, ^ being the chemical potential. 
The bandwidth is 2Ao = 2 max|^r(fc||)|- ffi and 172 are 
the backward and forward scattering amplitudes, respec- 
tively. We assume the band fiUing to be incommensurate 
and neglect Umklapp processes. 



A. One-loop RG equations 



The 2PI two-particle vertex *pp,pir''''"'''(^») now car- 
ries left /right indices Vi . Momentum conservation implies 
Ti = 0. We define 











= t,s). 




- $(2) + - 






= ch,sp). 




- $(2)++ 
A" 


~~(a;0, 




= ch,sp), (106) 



where C, P, and L refer to the Cooper, Peierls, and Lan- 
dau channels according to the standard terminology used 



in the framework of the g-ology modelii^ In Fourier space, 
we neglect the frequency dependence and approximate 

*Jr(A:i,fc2,fc3,^4) ^ <S>^{kF,-kF,-kF,kF) = <i>^, 

$J^(fci,fc2,fc3,/c4) ~ $P(fcF,-fcF,-fci^,fcF) = $1^, 
$J;(fci,fc2,fc3,/c4) ^ ^]:{kF,kF,-kF,-kF) = ^]l. 

(107) 



By scaHng arguments, one can show that the dependence 
of the IPI vertex 7^'*-'(fci) on uim and \ki\\ \ — kp is irrele- 
vant in the RG sense. The validity of Eq. Ill()7|l . which as- 
sumes that the 2PI vertex ^'^^^ shares the same property, 
will be discussed in Sec. IIV CI One-loop fiow equations 
for $^ and are deduced from Eqs. Ij54|l by includ- 
ing the r index, i.e. Xi (xi,ri). Retaining only leading 
logarithmic divergent loops (Parquet approximation) 
the contribution to the self-energy vanishes and the RG 
equations for the two-particle vertex read (see Fig. IHJ 



^ch = -i?pp($rh1>ch 



3$P $L \ 

'-^ sp spy 



-Br.r.(^K<^^ + <i>^ ^\ - 2$^ ) 

pp V ch^sp ^ ^sp^ch ^^sp^sp/J 



-2Sph 



ch 



3$ 



PP I 'J'ch 



3$^ 



*sp = -Sphl 4<i>t" 



c^c 



cn sp ' sp 



-5pp(*cWp-'i>s'p 



cn sp ' 



sp 



,(108) 



where 



^pi' = -^J2^G+{k)G^(k ~ Q) + [G ^ G)], 
^pp = ■^T.^G+ik)G^i-k) + (G ^ G)] (109) 
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come from the ph and pp loops, and Q = (2fcF,0). The 
dot denotes a derivation with respect to I = ln(Ao/A). 
To evaluate Sph and Spp, we use a sharp infrared cutoff, 



aw^cw^-^'l^-WI-f, (110) 



iuJn - S.r{k\\) 



This gives 



Spp = -5ph = ^tanh(4). (Ill) 

The initial values of the 2PI vertex is the bare vertex 
defined by Eq. 1110511 . 

*?|ao = -51 +32, $?|Ao = <?l+5'2, 
^chlAo =-f- +.92, $^p|Ao = -f. (112) 



B. Response functions and IPI vertices 

The response functions in the Cooper and Peierls chan- 
nels are defined by 

^ k,k' 



Xch,sp 



If E k-Q,k'-Q, fc').(113) 



BL ^ 

^ k,k' 



From the Bethe-Salpeter equations satisfied by xts ^ 
and xj)h~sn^ (Eqs. II59II with the index included), we 



obtain 



n 



pp 



" i + npp<i>c' 



-2n. 



ph 



1 - 2nph$p ' 



(114) 



where 



npp = ^E^+(^)G-(-^) 

= -Hpp. (115) 

Similarly, for the IPI two-particle vertex 7^'*) [Eq. l(fiT|l ]. 
we find 



7^ 



7.^ 



i + npp4>c' 



A' 



1 - 2nph$p ^ 



7^ = 



(116) 



The equality between 7^ and <I>J; is due to the absence of 
logarithmic divergent loops in the Landau channel and 
holds at the Parquet level. 

In Sec. m we have shown very generally that the RG 
equation for 7'-'*^ derived within the IPI scheme follows 
from the RG equation satisfied by the 2PI vertex <i>^^'' and 
the Bethe-Salpeter equation relating 7^'') and . In ap- 
pendix El as a means to check the validity of Eqs. IllflSII . 
we recover the RG equation satisfied by 7*^''^ directly from 
Eqs. 111(18111611. 



C. Discussion 

Figure El shows <i>^^^, x and 7^''^ versus I = ln(Ao/A) 
at zero temperature for gi = §2 = 0.2. gi — gi/nvp and 
52 = 32/7rwF are dimensionless coupling constants. Solid 
lines show the results obtained within the 2PI scheme by 
solving the fiow equations ljl()8|l and using Eqs. (jll4lll6|l . 
Dashed lines correspond to results obtained within the 
IPI scheme, where 7^^^ and x are directly obtained from 
RG equations, the 2PI vertex being then deduced 
from the relations II116|I . The agreement between the two 
schemes is excellent at high energies (small I), but deteri- 
orates at lower energies where, in at least one correlation 
channel, the 2PI vertex <i>'-^'' becomes of order one and 
eventually diverges. This deficiency is not important in 
the IPI scheme, as <i>*^^-' is usually not considered in this 
scheme, but it shows that the one-loop approximation 
breaks down in the 2PI scheme at low energy. 

The divergence of the 2PI vertex in the IPI 

scheme, which correctly predicts all other physical quan- 
tities, suggests that the relations II116|I between IPI 
and 2PI vertices may not be quite correct. Inverting 
Eqs. H116|l and considering the zero temperature limit, 
we obtain $^(0 = 7^(Z)/(1 - ^7^*^/2) in the Cooper chan- 
nel and ^^{l) — 7^(0/(1 ~ ^7^) in the Peierls channel. 
$ = ^/ttvf and 7 — 'f/nvF are dimensionless vertices. 
We conclude that the 2PI vertex will diverge at a 
finite energy scale Ac = Aoe~''= whenever 7*^'*-' is positive 
and finite in the low-energy limit {I 00). For repulsive 
interactions, 7^(0 ^ 52 — 31/2 satisfies this condition 
when 32 > .91 /2, which leads to the divergence of the 2PI 
vertex $^ (see Fig.|2l. The unphysical divergence of the 
2PI vertex <i>(^) comes from the assumption that the 2PI 
vertex is momentum independent [Eq. Hlfl7|l ]. which re- 
sults in an artificial decoupling of energy scales between 
the 2PI vertex and the reducible pp (or ph) propagator 
Hpp.ph- The latter involves all energy scales between A 
and Aq whereas the momentum-independent is an 
effective 2PI vertex at the energy scale A. The Bethe- 
Salpeter equations Ijll4|l and II116|I are therefore expected 
to be less and less reliable as the infrared cutoff A de- 
creases. 

These difficulties can be partially overcome by modify- 
ing the Bethe-Salpeter equation relating ^'^^^ to 7^'*-' and 
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FIG. 9: 2PI vertices (a), response functions x (b) and 
IPX vertices 7'''' (c) versus I — ln(Ao/A) for pi = 32 = 0.2 
and T — 0. Units are chosen such that ttvf — 1. The 
correlation channel is indicated in the upper left or right 
corner of each graph, following the notation of the text: 
(C,t) = (Cooper,triplet), etc. The solid (dashed) lines show 
the results obtained from the 2PI (IPI) RG scheme. In the 
IPI RG scheme, the 2PI vertex (dashed lines in panel 
(a)) is deduced from 7''''' using the relations Illl6|l . 
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FIG. 10: Sam e as Fig. HI but using Eqs. 111811 instead of 
Eqs. <114l11(it . 



7^(A) 



X^(A) 



<f,^(A) I " |tanh ( /3i ) l'^(A)7;r(A) 



^^(A) - I " |tanh (^/3f ) $^(07,^ (A), 



1 



npp(A) --J^ ^tanh j 4>^(A)x;r(A) 

- npp(A) - |tanh $^(e)x^(A), 

(117) 

and similar equations in the ph channel. The 2PI vertex 
^'^^^ is now taken at the same energy scale than the pp 
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or ph propagator. We thus obtain 



7^^ 



2 



C _ ^pp 



1 



IT/C 

2 



-2n 



ph 



where 

/a T 
satisfies the fiow equation 



tanh ( ] 



A 



A=? 



vi.^^P = tanh U^^P. 



(118) 



(119) 



(120) 



Fig. [mi shows $(2^ X and 7'^'*) obtained from Eqs. Ijll8|l 
instead of Eqs. II114I116II . The parameters are the same 
as in Fig. |2| We see that the agreement between the 
IPI and 2PI schemes, in particular for the susceptibil- 
ities, is significantly better. It is instructive to con- 
sider the case gi — where the IPI vertices 7^ and 7^ 
are fixed points of the fiow equations at any order in a 
loop expansion. From Eqs. Ijll8ll2()|l . one then obtains 
<i>c = a>C(? = 0)exp(i7^0 and = = Q)eMlll) 
at zero temperature. Depending on the sign of 7^, the 
2PI vertex will either vanish or diverge exponentially 
in the limit I ^ oo. The unphysical divergence of the 
2PI vertex obtained earlier at a finite energy scale is 
now replaced by an exponential divergence at zero en- 
ergy (A = Age"' — > 0). This divergence follows from 
the coexistence of logarithmically divergent (reducible) 
pp and ph loops and weak IPI vertices, which is made 
possible in ID by the strong interferences between var- 
ious correlation channels. The 2PI fiow equations give 
the correct qualitative behavior of '^'^'^\ but the one-loop 
approximation clearly breaks down at low energy when 
<i>*^2) becomes of order one. 

It should be noticed that if, within the IPI scheme, the 
susceptibilities were calculated from 



n 



pp 

2nph 



-^ppT;^ -^pp, 

4nph7?^nph, 



(121) 



with momentum independent vertices 7*^^-' obtained from 
RG equations, similar difficulties would arise and the 
correct asymptotic low-energy behavior would not be 
reached. The correct result is obtained by also deriv- 
ing RG equations for the susceptibilities; in this way ir- 
reducible and reducible pp and ph loops are considered 
on equal footing. These RG equations are obtained by 
introducing bosonic external sources that couple to or- 
der parameter fields, in addition to the fermionic sources 
that are used to obtain the generating functional of IPI 
vertices by a Legendre transformationiS2i2£ 

To quantify the failure of the 2PI scheme, within the 
one-loop approximation, to access the low-energy limit 



of ID systems, let us consider the case of the quasi-lD 
organic conductors of the Bechgaard salt family. In these 
systems, the bandwidth 2Ao ~ 4i|| ~ 12000 K is much 
larger than the kinetic interchain coupling ij_ ~ 300 K. 
For gi — g2 — 0.2, the 2PI fiow equations break down 
for Z ~ 5 (Fig. nrHl , which corresponds to an energy or 
temperature scale Age"' ~ 40 K. For 51 = 52 = 0.4, 
we find that the fiow equations remain valid down to 
^ Aqc"'^ ~ 300 K. In both cases, these energy scales are 
of the same order of magnitude or smaller than the tem- 
perature Tx ^ ^ 300 K at which a crossover to a 2D 
regime takes place. Preliminary calculations in quasi-lD 
syste ms ind icate that the interaction strength studied in 
Refs. IL'tIiH . namely 32 = 2^1 = 0.64 is accessible within 
the 2PI schemei^ We thus conclude that the 2PI scheme 
can be used for realistic quasi-lD systems like the Bech- 
gaard salts. 



V. SUMMARY AND CONCLUSION 

We have discussed the implementation of a Wilsonian 
momentum-shell RG approach within the 2PI formalism 
introduced in the 60s by Luttinger, Ward, Baym and 
others jiiSiSii The 2PI RG scheme yields an infinite hier- 
archy of fiow equations satisfied by the 2PI vertices 
The susceptibilities are obtained from the Bethe-Salpeter 
equation that relates them to the 2PI two-particle vertex 
<i>*^2^. In the normal phase, one has schematically 

App — APP 2 PP PPXPP 
\ 

(RG equations) 

K 

Xph = Xph + 2nph$phXph 

(similar equations could be written for the IPI vertices 
7pp^ and 7p^^ ) where ph are the susceptibilities with- 
out vertex correction and Hpp/ph the pp or ph pair prop- 
agator. The arrows indicate the coupling between the 
pp and ph channels which is taken care of by the RG 

f2) 

equations satisfied by the 2PI two-particle vertices $pp 

f2) 

and $pi/. The infinite hierarchy of RG fiow equations 

satisfied by the 2PI vertices should be truncated at 
some order in a loop expansion. The simplest non trivial 
truncation, the one-loop approximation, was discussed in 
detail in Sec.^ As any approximation of the 2PI vertex 
<i>(2\ it leads to a violation of the crossing symmetries of 
the two-particle Green function VF^^^ and the IPI ver- 
tex 7''*'. Besides, the one-loop approximation is not a 
<i>-derivable approximation, and a detailed study of con- 
servation laws and Ward identities remains to be done. 

We have shown in Sec. lIVI that ID conductors are char- 
acterized by an exponentially divergent 2PI vertex in the 
zero-energy limit, a consequence of the strong interfer- 
ences between correlation channels in ID. This leads to a 
breakdown of the one-loop approximation, which is there- 



19 



fore unable to access the asymptotic low-energy behavior 
of Luttinger liquids. Nevertheless, we have argued that 
the 2PI scheme can be used in quasi-lD systems like the 
organic conductors of the Bechgaard family where a di- 
mensional crossover always drives the system towards a 
2D or 3D behavior at low energy. 

On the other hand, the 2PI RG scheme leads to a par- 
ticularly simple description of single-channel (i.e. mean- 
field) theories. In mean-field theories, all the A depen- 
dence comes from the two-particle-reducible part of the 
susceptibilities (and the IPI vertex 7'^'), while the 2PI 
vertex is invariant under the RG transformation: 
$(2) = 0. The one-loop fiow equations then reduce to 
a single equation that determines the (anomalous) self- 
energy and reproduces the usual mean-field gap equation 
(Seeing. 

The possibility to continue the RG fiow into broken- 
symmetry phases is an essential feature of the 2PI RG 
scheme and is due to the fact that the 2PI two-particle 
vertex contrary to its IPI counterpart, is not sin- 

gular at the phase transition. This property, which is 
obvious in a mean-field theory, deserves some discussion 
in more general cases. The interchannel coupling, even 
weak, is likely to induce singularities in the 2PI vertex 
in particular in the vicinity of a phase transition. 
However, these singularities can be controlled by a proper 
parameterization of <i>^2)_ 'Jq ggg ^j^jg^ jg^ consider 
a quasi-lD conductor near a spin-density-wave (SDW) 
instability. The spin susceptibility Xsp(2fcF,7r) diverges 
at the SDW transition but, as in single-channel (mean- 
field) theories, the irreducible vertex <i>fp should remain 
finite (see Sec. IIVI for the definition of <I>^ and in a 
(quasi-) ID system). In the vicinity of the phase transi- 
tion, nearly divergent spin fiuctuations strongly affect the 
2PI vertex $^ in the Cooper channel. For the purpose 
of our discussion, let us assume that is proportional 
to the spin susceptibility, 

^M=t,s(fc-L>9-L - kj_,-k'j_,k'j_ - qj_) oc Xsp{2kF,kj_ +fc^), 

(122) 

where we retain the k± dependence of the vertexjifi The 
combination of k± arguments in Eq. II122II is the one 
that appears in the Bethe-Salpeter equation determin- 
ing the superconducting susceptibility. Expanding the 
even function Xspi'^kp, q±) in Fourier series, we obtain 

^^{k±,q± - k_L,-k'j_, fcl - q_L) 

00 

= ^a(") cosKfc^ + fcl)] 

n=0 

00 

= ^ a^"^ [cos{nk±) cos{nk'j_) — sin(nfc_L) sin(nfc^)], 

n=0 

(123) 

where a^""^^^a^"-' < and la^"-*! is a decreasing function 

of n. The condition la^T^I = \a^°''\, i.e. a|r^ = (-1)"4°\ 
would lead to a diverging susceptibility Xspi'2kF,q±) cx 



d{q± — tt), while ja^""''"'^'' | < |a|,"''| gives a broadened peak 
at q± = TT. The proximity of the SDW transition man- 
ifests itself by a larger and larger number of a^"'' coeffi- 
cients with a significant amplitude. All these coefficients 
should however remain bounded, with 1± (a^"''/2)npp > 
0, for the system to be stable against a superconducting 
instability (Hpp is the pp propagator defined in (II 1511 1. 
By parameterizing the 2PI vertex by means of the 

(n) 

all coefficients, we avoid any complication due to the 
diverging spin susceptibility. In practice, only a finite 
number of coefficients need to be retained. The coeffi- 
cients with a large value of n, which correspond to pairing 
between fermions n chain apart, do not play an impor- 
tant role; they cannot drive a pairing instability - an 
instability with a small value of n (n = 0, 1, 2, • • • ) will 
always occur firstSS - and their infiuence on the ph chan- 
nel is expected to be negligible. This approximation will 
lead to a non essential smearing of the singularity of 
at k± + k'j^ = TT. The argument given here for quasi- 
lD systems can be made more general. One can expand 
the 2PI vertex on the eigenmodes of the Bethe-Salpeter 
equation, retaining only a finite number of eigenmodes. 
Each coefficient in this expansion is bounded by a critical 
value at which a phase transition occurs. 

In Sec. Ill Fl we have proposed a method to derive the 
Ginzburg-Landau expansion of the thermodynamic po- 
tential in the vicinity of the phase transition on the basis 
of the fiow equations in the normal phase. This is the sim- 
plest way to access phases with long-range order as it does 
not require to solve the full RG equations in the presence 
of symmetry breaking. The knowledge of the Ginzburg- 
Landau expansion of the thermodynamic potential allows 
one to study the possibility of phase coexistence below 
the transition temperature. This is a particularly impor- 
tant issue in quasi-lD organic conductors where super- 
conductivity and antiferromagnetismj^Siffli as well as spin- 
and charge-density waveS(Si*£^ coexist in some regions of 
the pressure-temperature phase diagram. Recent IPI RG 
calculations have indicated that antiferromagnetism and 
superconductivity, as well as charge- and spin-density- 
wave phases, lie nearby in the phase diagram of quasi- ID 
conductors piS'i^'iS'ifi The 2PI RG scheme would allow 
to determine whether they actually coexist in the low- 
temperature phase as observed experimentallyi^SiSSiSiiS^ 

Finally, we note that the 2PI RG scheme enables a di- 
rect connection to the phenomenological Landau Fermi 
liquid theory when the metallic state remains stable down 
to low temperature. It has been shown in Ref.lfi^that the 
functional T[G], or rather its variation ST — r[G] — T[G], 
can be written as a functional (5r[(Sn] of the Wigner distri- 
bution function n — {nk(T(r, t)} which is essentially de- 
termined by the Landau parameters. The latter are given 
by the forward-scattering Hmit of the 2PI vertex <i>(2). 
{ftn = n — n denotes the deviation from the equilibrium 
value n.) 6T[Sn] determines both static and dynamic 
properties of the Fermi liquid. In particular, it yields the 
quantum Boltzmann equation satisfied by nkcr(r, t). In 
the static case, the Wigner distribution {nk(T(r,T)} re- 
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duces to the quasi-particle distribution function {rikcr}, 
and 6r[Sn] = SE[6n] - fiSN - TSS[Sn] where 6N[Sn] 
and <55'[(Sn] are the quasi-particle number and entropy 
variations induced by a change Sn in the quasi-particle 
distribution. The functional SE[Sn], which gives the cor- 
responding energy variation, is the basis of Landau's phe- 
nomenological Fermi liquid theory. 
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where we have used Bpp = — Bph. 7 , 7^ and are 
defined similarly to <i> , (f>^ and <1>^ [Eqs. 1110611 ]. Note 
that 7('')ri7('*) vanishes in the Landau channel as it does 
not produce any logarithmic term. Eqs. (|A1|I can be 
simplified by using the crossing symmetries 
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APPENDIX A 



In this appendix, we show that Eqs. Ijl08|l correctly 
reproduce the one-loop equations for the IPI vertex 
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in Luttinger Hquids. The flow equation of 7 



(4) 



is given by Eq. ll39)l . To order (7^**^)^, it reduces to 
7(4) ^ _7(4)n^(4) _^^(2) ^(2) should be evaluated 

to 0[(7(4))2]. This leads to 
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From Eqs. |IA1IA3|I . we finally deduce 
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with the initial value 7 
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(A4) 



$(^)|ao given by 



Eqs. II112|I . Alternatively, using the parameterization 
H112|l . one can rewrite Eqs. I)A4|I as two flow equations 
for 51 and 52, 



51 = -2Bpp5i, 

252-51 = 0. 



(A5) 



Eqs. IIA5|I are the usual one-loop RG equations for a ID 
systempi^ 
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